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Abstract 

We perform a calculation of the first and second order infinitesimal variations, with respect to 
energy, of the Boltzmann entropy of constant energy hypersurfaces of a system with a finite 
number of degrees of freedom. We comment on the stability interpretation of the second vari- 
ation in this framework. 
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1. Introduction 



In theories describing systems with many degrees of freedom, field theories being an exam- 
ple, it is of great interest to be able to derive thermodynamic quantities from the microscopic 
dynamics. The microscopic dynamics of a non-dissipative system is encoded in its Hamiltonian 
description. In a statistical description of such a Hamiltonian system, we trade the practically 
intractable symplcctic evolution on its phase space with a "reasonable" probability measure 
that describes some of the characteristics of such an evolution. In systems in thermodynamic 
equilibrium the use of the microcanonical and the canonical distributions has proved exceed- 
ingly successful during the last century and these arc the distributions with respect to which 
we will be calculating the statistical averages in this paper. 

The Boltzmann entropy has proved to be one of the most useful thermodynamic potentials. 
The Boltzmann entropy is proportional to the area of the total energy E hyper-surfaces Ai e 
on which the configuration space M of the Hamiltonian system can be foliated. Because of this 
very direct geometric interpretation, we have chosen to analyze variations of the Boltzmann 
entropy. Naturally, the different response functions of such a thermodynamic system can be 
derived in terms of appropriate variations of the Boltzmann entropy. The most geometrically 
transparent, and at the same time physically relevant, of such variations are the ones with 
respect to the total energy E of each hypersurface, to which we focus. 

An outline of this paper is as follows: in Sections 2 and 3 respectively, we calculate the first 
and second order change of the Boltzmann entropy under infinitesimal variations of the energy 
E (diffeomorphisms) of M.e- We also provide a physical interpretation of these results in 
terms of models without kinetic terms (e.g. lattice models). In Section 4 we present some 
conclusions and the comment on the relation of this paper with similar works. 

2. First order variation of the entropy 

Assume, following the ideas of Krylov [1] that A/", the configuration space of the system 
under study, is an n-dimensional Riemannian manifold with metric g. The evolution of such 
a system is described by the geodesic fiow of g. Consider two diffeomorphic, infinitesimally 
close, hyper-surfaces M.e and A^e+^b of H with corresponding energies E and E + SE 
respectively. Since the system is autonomous, its evolution can be described by restricting our 
attention to {A4 e^q)- When the system is coupled to a heat reservoir, because of the existence 
of energy fluctuations 5E the system may flnd itself in an "adjacent" hypersurface M.e+5e- 
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If 5E <^ E and the system is not close to a phase transition, then it is reasonable to expect 
that A4e and A4e+se are diffeomorphic. 



The thermodynamic potential with the most straightforward geometric interpretation is the 
Boltzmann entropy. The Boltzmann entropy describes accurately the behavior of a system, as 
long as the system evolution is ergodic [2] on its configuration space Af. Since it is difficult 
to prove, in practice, the ergodicity of a system of physical significance starting from first 
principles, we proceed by assuming that the geodesic flow is ergodic [2] on both Me and 
M.E+SE when they are diffeomorphic to each other. Let f ■ M.e ^ M-e+se be such a 
diffeomorphism, and X e TM\m be the vector field generating /, when restricted on M.. 
X can be written a^s X — and let £x denote the Lie derivative along X. We assume 

E 

that the metric on AIb is the induced metric g from g [3]. Let V and V indicate the 
Levi-Civita connections on H and Ai compatible with g and g respectively [3] . Then 



(1) 



All Latin indices take values 1 through m in the sequel. Although M.e is obviously a 
codimension 1 sub-manifold of A/", during most of the calculation we will be more general, 
namely we will assume that the dimension of M.e is m < n without necessarily m = n — 1. 

As is well known, the Boltzmann entropy of AAe is given by 



Se — ks In volAA E 



(2) 



where denotes the Boltzmann constant and 



volM E = J Vdet g dTx 

Me 



Then 



5S := Se+5e - Se ^ ks In 



voIMe+se 



voIMe 

By keeping terms up to second order in 5E and upon expanding the logarithm we find 



(3) 



(4) 



5S 



d{volMi 



volM E dE 



6E + 

^ 2voIMe 



(P{voIMe) 



dE^ 



1 fd{volME) 
voIMe \ dE 



E 



{5Ef 
(5) 



which can be re-expressed [3] as 
5S = -4^{^xvoIMe)SE + 



voIMe 



2 voIMe 



Zx^^xvoIMe ^ {^xvoIMe? 

voIMe 



{6Ef (6) 
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Let {ci}, i — 1, . . . ,m be an orthonormal basis of TA4 with respect to g. Let V, V 

denote the Levi-Civita connections compatible with g and g on. H and Ai respectively 
[3]. Then (1),(3) give 

SIxvoIMe^ J £x^detg{ei,ej) d^x (7) 

Me 

For any positive-definite symmetric matrix A, we use the identity del A — exp(Trln74), where 
Tr denotes the trace of A, and we get 

2.xvoIMe^\ j ^det^(ep,e,) S.xTrhi[g{ei,ej)\ dTx (8) 

Me 

We observe that the Lie derivative £x and the trace Tr operations commute, and by using 
(1) 

£.xvoIMe = \ j y^^M^^ TrZxH9{ei,e,)](rx (9) 

Me 

which eventually gives 

2xvoIMe^^ j ^Jdetg{ep,eg) [g{ei,ej)]'^£x[g{ei,ej)] dTx (10) 

Me 

where repeated indices are summed. Since g{ei,ej) e R are functions on Me [3] 

S^x[g{(^i,ej)] = Vx[g{(^i,ej)] (11) 

Then 

Vx[^(ei,ej)] = {Vx9)iei,ej) +g{Vxei,ej) +g{ei,Vxej) (12) 

Since V is the Levi-Civita connection compatible with g, then [3] (Vx^)(ei, Cj) = and 
its torsion is zero, namely 

Vxei - Ve,X = [X, ei] (13) 

where [X, e^] denotes the Lie bracket between the vector fields X and e^. The metric 
g gives rise to the orthogonal decomposition X — + X-^ into a tangential component 
X'^ e TM.E and into a normal component X-^ e NM.E, where NM.E denotes the normal 
bundle of M.e- Then (12), (13) and linearity imply 

^xWi, ej)] = g{Ve,X^, ej) + ^([X^, e^], e,) + ^(V^xe^, e,) + {i ^ j) (14) 

where {i <-> j) indicates similar terms with i and j interchanged. The first term of the 
above sum is the second fundamental tensor [3] lx^{e-i-i ej). By using (13) again, we get 

Vx[^(ei, ej)] = lx^{ei, ej) + g{[X, a], ej) + gCVe^X^, ej) + {i ^ j) (15) 
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Substituting (15) into (10) we find 

^xvoIMe= j ^det^(ep, e,) [/^^(ei, e^) + ^([X, e^], e^) + ^(Ve,^^, e^)} (16) 

Me 

The determinant under the radical is equal to one, since {e,} is an orthonormal basis with 
respect to g. The third term is, by definition, the divergence of . Then by Green's 
theorem, we find 

y"^(Ve,xT,e,)= j g{X'^^y)dii (17) 

M.E OMe 

where v e TdM.E represents the outward unit normal on the boundary SAAe and is 
the induced Riemannian measure on OMe- If X is perpendicular to M.e, i-e. if X~^ is 
zero, or if Me is closed, then this term is trivially zero. Let Cj(s) denote the integral curve 
of ej, i.e. ^^^7^ = 62 and let Ci{s,E) be the one parameter variation of Ci{s) along X, i.e. 
Ci{s, 0) = Ci{s) and ^ x. Then 



d_ 
ds 



Ci{s,E)^eiX (18) 



so the second term of (16) is zero. If the above conditions hold, then (16) simplifies to 

£.xvoIMe^ j ydet^(e^^ lx-L{ei,ei) d'^x (19) 



Mi 



Let (3 — l/ksT, as usual. If the system under study has a constant extensive variable V, 
e.g. volume, and constant "particle number" J\f, then 



dS 
dE 



= ^ (20) 



and using (6), (20) we obtain 



P = — -r , (21) 

j y/det g{ep, Cg) d"^x 

Me 

which can be interpreted as the average of the mean curvature Tr Ix^ with respect to the 
micro-canonical measure 

and (21) can be rewritten as 

/3 = ( Tr /xx ) (23) 
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where the average ( ) is taken over Me- Since /3 > then ( Tr Ix-l ) > 0. An example 
where this condition is satisfied is when J\f — M.^, with A4 being diffeomorphic to the sphere 
S"~^ and isometrically embedded in R". Then Ix-l > everywhere on Al, a fact which 
clearly guarantees the positivity of (23). For this example, and because is non-compact, 
we assume that either — 0, or all the functions on A4e have compact support. Generally, 
however, we cannot exclude the possibility ( Tr Ix-l ) < 0. In such case, (23) loses its direct 
physical interpretation. One reason why ( Tr Ix-l ) < 0, can be traced to the lack of ergod- 
icity of the geodesic fiow on Me, which was assumed at the outset. Without such ergodic 
behavior, the expression for the Boltzmann entropy (2) is reduced to just a formal definition 
devoid of any physical meaning. This lack of physical meaning is subsequently inherited to 
thermodynamic relations like (20), where T can no longer be identified with the physical 
quantity "temperature". An alternative interpretation of (23), is as a constraint equation on 
the possible choice of a metric g describing the evolution of the system on Af. In such 
an interpretation, a metric g resulting in lx± < is not acceptable, on physical grounds. 
Therefore, either the metric g used for the description of the system should be modified, or 
in extreme cases, one should take the more radical step of discarding the model altogether. 

If the model under consideration, however, describes a system (lattice models are, frequently, 
such examples) for which there is an upper bound in the possible energy, then negative tem- 
peratures are theoretically allowed, in the definition of the partition function of a canonical 
treatment. Such models should not, obviously, contain any kinetic energy terms and this is 
refiected on the choice of ^ on jV describing their evolution. Systems being described by 
such models have been experimentally observed [4] to be out of equilibrium, a fact which puts 
in question, the suitability of the Boltzmann entropy in describing them, especially when they 
are far from equilibrium. For such models there is no constraint on the sign of (23), but to 
acquire a physical meaning, T should be interpreted appropriately. 



3. Second order variation of the entropy 

For systems in which both positive and negative temperatures have physical meaning [4], 
a limiting case occurs when the microcanonical average (23) of the mean curvature is zero, 
( /x-L(ej, Cj) ) =0, which amounts to (3 — or T being infinite. This requirement is trivially 
fulfilled [4] when Me is a totally geodesic submanifold of jV, i.e. when lx-L{ei,ej) = or 
when Me is a minimal submanifold of J\f, i.e. when the mean curvature lx-i-{ei,ei) ~ 0. 
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In such cases the Boltzmann entropy remains invariant under the action of X and (6) gives 

5S = ^ ^f. . (S^xS^xvoIMe) {SEf (24) 
To perform this calculation we start by Lie-differentiating (16) 

£x^xvoIMe = j !^S.x^Jdetg{ep,eg)^ l^lx±{ei,ei)+g{[X,ei],ei)+g(Ve^X^,ei)^ dTx^ 



E 



j {£x{/x4e^,e,)} + £x{^([X,e,],e,)} + £x{^(Ve,XT,e,)}} cTx 

Me 

The first term of the right hand side is given by (16). Taking into account (18), the definition 
of the second fundamental form, and that X = X"^ + X-^, we find 

ZxSIxvoIMe = j {^(Ve,X^, e,) + ^(Ve,X^, e,)}' ^/dit^ (Tx + 

Me 

j { Vx{^(Ve,X, e,)} + Vx{^([X, e,], e,)}} ^d^ cTx 



Ml 



which gives, after using the torsion-free condition (13) with (18) and performing the covariant 
differentiations 

S^x^xvoIMe^ j {{^(Ve,X,ei)f + ^(Vxei,Vxei)+^(VxVxei,ei)} ^/d^g (Tx (25) 

Me 

By using that V is Levi-Civita with respect to g, this can also be written as 

ZxS^xvoIMe^ j {{^(Ve,X,e,)F + Vx{^(Vxei,eO}}A/dit^d"^x (26) 

Me 

If M.E is a minimal submanifold of then (26) reduces to 

Zx^xvoIMe^ j [{g{VeX.^i)Y + ^x{g{Vxei,ei)}]./^g(rx (27) 

Me 

A further simplification occurs when X is everywhere normal to Ale, i.e. when X'^ — 0. 
Then 

£.x2.xvoIMe^ j Vx{g{Vxei,ei)} y/d^ d'^x (28) 

Me 

Substitution of (28) into (24) gives 

= ^ ^ff/ / Vx{^(Vxe„ e,)} d^x (29) 



M 



e 
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which can be re-expressed, by using (22), as the microcanonical mean 



^^^M^ (Vx{^(Vxe„e,)}) (30) 



or, equivalently, as 



^ = Y ^ ^(Vxe„ Vxe.) + ^(Vx V^e,, e,) ) (31) 
Since g has positive signature, the first term of the right hand side of (31) is positive or zero. 
For the same reason, the operator in the second term is eUiptic, thus it eventually has positive 
eigenvalues. Then each side of (31) can either be positive or negative, in general. Because of 
(20), and since the heat capacity Cy under the constant extensive variable V is 

dE 



(32) 



(30) can be re-written as 



^ = -^ (Vx{^(Vxe.,e,)}) (33) 

During the second order variation (3 = which, according to (33), implies that cither = 
or ( Vx{5(Vxei, e,)} ) = 0. In order to avoid extending the expansion (6) to cubic and higher 
order terms in SE, we consider only the former option. The result of (33) has a physical 
interpretation as long as 

is finite. If ( Vx{^(Vxej, Ci)} ) < then Cy > which is the standard stabihty criterion. 
On the other hand, if ( V x{g(^ x^i, Cj)} ) > 0, then Cy < which indicates that the system 
is unstable. We can, therefore, interpret ( V x{g(V x^i, Cj)} ) as a quantitative measure of the 
instability of a system with an upper bound on its energy. 



4. Discussion and conclusions 

Some of the above results are standard in the theory of minimal submanifolds [5], [6], [7] and 
the theory of harmonic maps [3]. In the second variation, we deviated considerably from the 
established practice [6], [7] which results in an inner product of X with an elliptic operator 
(Jacobi operator) expressed in terms of the Laplacian and of the Riemann tensor of the normal 
bundle NAi e acting on X. We did so because we did not need the aforementioned geometric 
result in order to obtain a physical interpretation of the second order variation of the entropy. 
Evidently our result can be recast in the form provided by [5], [6], [7] upon integration by parts 
and by using (23) with j3 = 0. 
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It may also be worth noticing the similarity of the present results to the ones of [8], [9]. 
In these papers the author relies mostly on measure-theoretical arguments to draw his conclu- 
sions. The use of a Euclidean metric on the phase space is very minimal [8] to none [9] . Clearly, 
measure-theoretical arguments [2], [9], [10], [11] are applicable to a much wider variety of systems 
than the mechanical Hamiltonian ones that we use here. Dissipative systems [10], [11] are an 
important class of systems that the Riemannian approach, as used in the present paper, cannot 
describe. On the other hand, the Riemannian approach may shed some light into aspects of 
Hamiltonian systems as, for instance, the relation between Gaussian curvature and dynamical 
temperature [8] which may not be so clear, or accessible, if one uses purely ergodic arguments. 
Such a relation has been pointed out by the author of [8], who curious as he was about it, made 
no attempt to trace its origins. Whether such a relation actually exists and can be elucidated 
by using Riemannian methods can be a topic of future research. 
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